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It is known that elementary bosons condense in a unique state, not so much because 
this state has the lowest free particle energy but because it costs a macroscopic 
amount of energy to put the particles into different states which can then interact 
through quantum particle exchanges. Since individual exchanges between the two 
fermions of a composite boson are ignored when composite particles are replaced by 
elementary bosons, it is of importance to reconsider the exchange energy argument 
for the stability of the Bose-Einstein condensate in the case of composite bosons. 
We do this here in the light of the new many-body theory which allows us to take 
exactly into account all possible exchanges between the fermionic components of the 
composite bosons. We confirm that the condensate of composite bosons is made of a 
unique state, this state being moreover pure: a coherent superposition of states close 
in energy is shown to be less favorable for both elementary and composite bosons. 

PACS numbers: 03.75.Hh, 67.10.Ba, 71.35.Lk, 67.85.Bc 



I. INTRODUCTION 

Bose-Einstein Condensation (BEC) is of high current interest due to recent observations 
of condensates made of atom^, molecule d^ * ^ * ^ * ^ !, and polaritond^, which result from the 
strong coupling of a photon and an exciton^. The experimental observation of a pure ex- 
citon condensatd^'^, however, remains a challenge. It has been recently argued^^ that the 
condensate in this system should appear in a dark state, since these are the lowest energy 
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states due to the weak valence- conduct ion repulsive processes which do not exist for excitons 
with spin S=(+2,-2). Consequently, such a dark condensate cannot be directly seen by op- 
tical emission investigated up to nov^, though it should be possible to deduce its presence 
indirectly. In addition, excitons are excited states; so that, in order to reach condensation, 
excitons with lifetime long compared to their thermalization time are needed. This has mo- 
tivated the development of coupled quantum well structured with electrons well separated 
from holes to increase their recombination time. 

The standard formulation of BEC considers a set of free elementary bosons. For particles 
with a center of mass momentum, a special role is then played by the zero momentum 
[k = 0) state. However, since the energy spectrum of such particles confined in a large 
volume is essentially continuous, it is physically hard to accept that this lowest energy state 
is favored over all the other nearby states just on account of an infinitesimally small kinetic 
energy difference. 

Actually, the essential characteristic of a Bose-Einstein condensate, that there is a macro- 
scopic number of particles in a single quantum state, arises not from an argument for an 
ideal gas but from interactions. This has been made in a very straightforward way by 
Noziered'^ By calculating the energy in the Born approximation of elementary bosons 
with repulsive interaction (necessary to avoid a density collapse) he has shown that to break 
up the condensate into two different states, we must pay a macroscopic exchange energy 
penalty which increases as A^^. As Nozieres said, "...it is the exchange interaction energy 
that makes condensation fragmentation costly. Genuine Bose-Einstein condensation is not 
an ideal gas effect: it implies interacting particles!" 

This nicely shows that exchange between undistinguishable quantum particles plays a 
fundamental, not peripheral, role in the essential physics of Bose-Einstein condensation. A 
similar role is actually played in the spin-spin J-coupling of ferromagnets, favoring particles 
to be in the same state: Without this coupling, there is no ferromagnetic phase transition. 
As shown below, even for just N = 2 bosons, the result is still valid; in other words, there 
already is an exchange energy penalty for two bosons to be in different states. As a general 
rule, we can say that Bose-Einstein condensation occurs when the exchange energy for two 
bosons becomes comparable to their thermal energy. 

As Nozieres' conclusion was obtained before the development of the composite boson 
many-body theory, it was de facto reached in the framework of elementary bosons. Yet, all 
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real condensates consist of composite bosons which are made of an even number of fermions. 
Due to the Pauh exclusion principle between fermions, additional exchange processes be- 
tween the fermions which make up the composite bosons must be considered in the overall 
exchange energy. These effects are expected to be important not only for excitons composed 
of light-mass electron and hole, but also for atoms, since the electron exchange energy is 
known to enter the effective interaction between two atoms. 

A number of work^^ ^ ' ^^ ' ^^ ' ^° ' ^^ ' ^^ ' ^^ ' ^"^ ' ^^ ^^have addressed interaction and condensation of 
composite bosons ("cobosons" in short) such as excitons. One approach is to stay completely 
in the fermion picture. Treating the correlated pairs in this approach can be cumbersome 
and can require heavy numerical methods. Another popular approach is the method of 
bosonization, in which the fermion pairs are treated as pure bosons but with an altered 
interaction taking into account the underlying Fermi statistics. As discussed in a series of 
recent papers presenting a new approach to composite boson theorj^^^l^, the bosonization 
method neglects certain exchange processes which are important even in low-order perturba- 
tion theory, and neglecting these can sometimes have dramatic consequences, for example, 
neglecting dominant terms in semiconductor optical nonlinearitieP^'^. In the present paper 
we will use the new composite boson theory, which has a convenient diagrammatic method 
which lends itself to analytical results. We reconsider the overall stability of Bose-Einstein 
condensation in the case of composite bosons using the framework of this new many-body 
theory, since exchange, which is crucial for the condensate stability, is not properly treated 
when the fermionic components of the particles are forgotten. Our aim in this paper is 
to show that, indeed, the exchange-energy stability argument still applies in the case of 
composite bosons, using the new tools that this composite boson many-body theory now 
offers. 

This many-body theory shows that two composite bosons interact through two conceptu- 
ally different scatteringd^ "interaction scatterings" for fermion interactions in the absence 
of fermion exchange, which have energylike quantities, and dimensionless "Pauli scatterings" 
for fermion exchanges in the absence of fermion interaction. These Pauli scatterings, by con- 
struction ignored when the composite particles are "bosonized", turn out to be crucial in 
the many-body physics of composite bosons: they, in particular, control all semiconductor 
optical nonlinearities induced by unabsorbed photons. From these 2x2 Pauli scatterings, 
we can construct any possible fermion exchange which exists between N composite bosons. 
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as necessary since the Pauli exclusion principle from which these exchanges originate, is 
A^-body by essence. These A^-body exchanges are nicely visualized through new diagrams, 
called "Shiva diagrams"'^, which not only allow one to see the subtle many-body physics 
taking place between these tricky objects but also to calculate it readily. 

This new formalism also shows that it is impossible to write an effective Hamiltonian for 
bosonized excitons which produces the correct scattering rates and the correct lifetime of 
exciton stated, even in the extreme dilute limit of = 2. A way to grasp the difficulty is 
to note that, by mapping composite bosons into an elementary boson subspace, we strongly 
reduce the degrees of freedom of the problem. This mathematically shows up through a 
change from 1/A^! to (1/A^!)^ in the prefactor of the closure relation for A^ elementary or 
composite boson^, making all sum rules irretrievably different. 

In this paper, we consider the possibility to have not only a condensate made of two 
different states but also a condensate made of a coherent superposition of states close in 
energy, for both elementary and composite bosons. Thus, we are going to consider the three 
states, 

i0o) = io), (1-1) 
i0i2)=5r^if io), (1-2) 

\<P) = {a'Bl, + a"Bir\Q), (1-3) 

where the B\ operators are creation operators for composite bosons, and we will compare the 
results with the ones obtained with elementary boson operators, written with a bar as B\. 
These operators having the standard commutation relation [Bm,B\] = 5mi- For composite 
bosons, the B\ operators obey the relation [Bm, B\] = 6mi — Dmi, where the "deviation-from- 
boson" operator Dmi will be discussed in more detail below. The third state, defined 



in (1-3) is a coherent superposition of composite bosons in two different eigenstates. The 
possibility to have such a coherent superposition, not considered by Nozieres but considered 
by Leggettf^, will also be addressed here in the case of elementary bosons. In the spirit 
of the Nozieres and Leggett arguments, we consider only the lowest order of perturbation 
theory, calculating the mean-field energy of the different states. This is certainly valid in 
the low-density limit, although as discussed by LeggetlP^, a rigorous demonstration of the 
stability of a condensate in the general case of interacting bosons at all densities does not 
yet exist. We confirm that the mean value of the Hamiltonian is minimum for |0o) and \4>o)] 



5 



i.e., for a condensate which is pure and not fragmented, whatever the repulsive scatterings 
between particles, and whether or not the bosons are elementary or composite. 

As the calculations for a large number of composite bosons are heavy and quite technical, 
especially for readers who are not yet familiar with this new many-body theory for composite 
bosons, we first perform the calculations for N = 2. Most of the important physics is 
usually seen just from examining two interacting particles. A nice rule of thumb — which 
will be explicitly confirmed — allows us to obtain the results for by replacing [as/LY in 
the results for A^ = 2 by 

r/ = N{aB/LY, (1-4) 

where is the composite boson wavefunction extension, L the sample size, and d the 
spatial dimension. This rule of thumb is less obvious for the fragmented state |0i2) because 
we then have two boson populations A^i and A^2- However, we physically expect to have the 
exchange terms between the composite bosons in states 1 and 2 to appear with a prefactor 
N1N2, while the ones involving exchange only in state 1 or only in state 2 should appear with 
prefactors Ni{Ni — l)/2 and A"2(A"2 — l)/2 for the number of ways to choose two particles 
among the A^i or A^2 populations. 

As with other many-body effects between composite bosons, the tricky part of the cal- 
culation always is to determine the density expansion of the relevant scalar products for 
A^ composite boson states. While it is always possible to get them through a brute force 
algebra making use of the commutation^ on which the composite boson many body theory 
is based, these expansions are nicely performed by using Shiva diagrams'^^ which visualize 
fermion exchanges in a transparent way: the density expansions are naturally associated 
with diagrams having an increasing number of composite boson lines. In the case of state 
1012), however, we must work a little harder, because we now have two large numbers which 
are relevant, A^i and A^2- In order to determine the extra energy due to the fracturing of the 
condensate into two different states, we must carefully distinguish between the exchanges 
which are internal to the A^i population, those which are internal to the A^2 population, 
and those which occur between these two populations. Only these last contribute to the 
exchange penalty for fragmenting the condensate into two different states. 

The paper is organized as follows: 

In Sections [TT] and III, we review the formalism for many-body effects with elementary 
and composite bosons, and settle the notations. An important part of this discussion is the 
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normalization of A^-particle states, which is trivial for elementary bosons but far from trivial 
for composite bosons. In particular, we present new results for the normalization of a state 
made of two large numbers of different cobosons. 



In Section IV, we calculate the mean value of an effective Hamiltonian describing inter- 
acting elementary bosons, in the case that of two elementary bosons in each of the three 
configurations |0o), \4>i2) and while in Section |v] we perform the same calculations 
for two composite bosons, using the exact Hamiltonian for interacting fermions. This in 
particular, allows us to identify the proper effective scatterings one has to take for diago- 
nal processes between bosonized particles in terms of the interaction and Pauli scatterings 
appearing in the composite boson many-body theory. 



In Section VI , we calculate the effective Hamiltonian mean value for elementary bosons 



in the three states of interest, namely pure, fragmented and coherent, while in Section [VII 
we address the same problem for composite bosons. 



In Section VIII, we discuss the results and conclude. 

In the appendix, we briefly reproduce Nozieres's original argument, for completeness, 
and also because the spirit of the present paper is the same: by calculating the Hamiltonian 
mean values in the states \4>o), \4>i2) and we by construction study the effect of inter- 
actions between bosons in the Born approximation, i.e., to first order in the interactions. 
States Bq^\0) and Bq^\0) are the system ground states in zero order for the elementary and 
composite bosons, respectively. For those states, the mean value of the Hamiltonian reduces 
to NEq in the low-density limit. By drawing our conclusion about the non-fragmentation 
of the condensate from the sign of the term linear in density, obtained within the Born 
approximation, we of course implicitly assume that higher-order terms in density will not 
modify the overall sign, which is likely in view of our past knowledge of many-body effects. 
This approach assumes zero temperature and low density. Actually, low density is implicitly 
assumed when considering excitons, which can convert to electron-hole plasma as density 
increases.'^ 
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FIG. 1: Scattering ^(^i)' defined in equation (2-2), between two elementary bosons. 



II. FORMALISM FOR ELEMENTARY BOSONS 

In this section, we briefly recall the formalism for elementary-boson many-body calcula- 
tions. The commutation relations for these bosons are 

[B^,Bi]=0, [E^,El] = 5„i. (2-1) 

The index i usually stands for a center-of-mass momentum Qi and a relative motion index 
Ui. For elementary bosons, it is possible to split the system Hamiltonian a.s H = Hq + V 
where the one-body part is Hq = J2 EiBjEi while the interaction can be written as 

V=lj:^ili)BlBlB,B,. (2-2) 

The energy-like prefactor C,{mi)j which describes the scattering from i to m and j to n, (see 
Fig. must be such that ^(^ i)* = ^(1 m) insure hermiticity. Momentum conservation 
requires that this scattering differs from zero only for + Qn = Qi + Qj- Nozieres' 
argumentl^l was made for a structureless scattering, i.e., for ^(J^l) = ^o^q^+q^ q^+q.- By 
noting that Bj^B}^ = B^B]^, it is possible to rewrite this potential \^ in a fully symmetrical 
form 

^ = 2 ^ ^'^'''^i^'^rnBiBjBi, (2-3) 



where imn;ij, equal to [^(;; ^) + ^(™ |)J /2, is now such that imn;ij = imn-ji = ^nm;ij = ilj-mn- 

In order to perform many-body calculations with elementary bosons in a convenient way, 
let us note that 

[Brn, 5^] = 5pf -1 + 5t[5„, 5pf -2 + ... = iV5_5f (2-4) 



This gives for the state |0o) defined in (1-1) 



S,-S,|0o) = N{N - l)5,„5,„5f -2|0). (2-5) 
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If instead of a single state, we now consider a coherent superposition of two states, 5^ 
a'Bl, + a" B\„, equation (2-4) implies 



= AT^.fitA^-^ (2-6) 

with 5i such that [Bi,B'^] = 6i, i.e., 

6i = a6io' + a"5io" (2-7) 

III. FORMALISM FOR COMPOSITE BOSONS 

A. Elementary scattering 

We now consider composite bosons made of one fermion a and one fermion (3, and we 
briefly review the formalism presented in Refs. |23and|2Hl The Hamiltonian of these fermions 
reads 

H = H^ + Hp + Vaa + Vpp + Vo^p. (3-1) 

For excitons or hydrogen atoms, Vaa^ Vpp and Vap are Coulomb potentials, while for the 
so-called cold Fermi gases, Vaa ~ ~ while Vap is short range. These interaction terms 
can be formally written as 

v^^-Y.vi!l,Y]<i.>i»A (3-2) 



a a 



'13 -13 



and similarly for V^a and Vpp with a 1/2 prefactor. at and fot are, for simplicity, chosen to 
be the creation operators for one fermion a or /? in an eigenstate of the system Hamiltonian, 
namely \ka) = ot |0) and \kp) = 6t |0), with {H — er )\k~f) = for 7 = a or /3. 

The creation operators of the one-pair eigenstates \i) = bI\0) of the system Hamiltonian, 
with [H — Ei)\i) = 0, can be written in terms of these free fermion operators as 



^1 = E«L4,(^/3^"I^)- (3-3) 
In the same way, the fermion pair creation operators can be written in terms of the -B|'s as 

{kpkali) is the wave function of the eigenstate i in momentum space if {Ha, Hp) only con- 



tain free energy contributions. These two equations (3-3) and (3-4) are quite fundamental 
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in the composite boson many-body theory as they allow us to "open" the composite bosons 
into their fermionic components to let the fermions interact through their exact interac- 
tion potentials and then to "close" these fermions back into composite bosons after their 
interactions. 

The i?|'s are composite boson operators. Indeed, their behavior is like bosons, since 

[Bra, 5,] = 0. (3-5) 

However, since these operators are such that 

[Bm, 4] = 5rm " I^™, (3-6) 

they differ from elementary bosons, the deviation from elementary boson statistics appearing 
in the operator Dmi- This operator is such that -Dmi|0) = 0, as can be seen by taking the 
above equation acting on the vacuum state |0). The "Pauh scatterings," which describe 
fermion exchanges between two cobosons in the absence of fermion interaction, are obtained 
from these deviation-from-boson operators Dmi-, through the commutation relation 

= {KJ^ + n)]Bl (3-7) 

where the second term means the same as the first but with m and n interchanged. A(^^), 
which is a dimensionless factor, describes the fermion exchange between cobosons («,j) 
shown in Fig. [2]^a). 

To get interaction scatterings which are energylike quantities, we introduce the "creation 
potential" ^ which arises from the commutation relatioEp3351 

[H,B\]=E,B\ + V^, (3-8) 

so that V^^|0) = for -BJ|0) being the H eigenstate with energy Ei . The interaction 
scatterings ^(^1) then appear through the commutation relatioEp^ 

[V^,B}]=Y.a'^l)BlBl (3-9) 

The prefactor ^{^i), shown in Fig. [2|b), describes the direct scattering of the "incoming" 
cobosons resulting from the interactions of their fermions, the "outgoing" cobosons 

{m,n) being made with the same fermion pairs. From 'C(mi)' possible to construct 
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FIG. 2: a) Pauli scattering A(5^^), as defined in equation (3-7), for fermion excliange between 
two composite bosons in the absence of fermion interaction. Note that this scattering can be seen 
as an exchange of fermions (3 or an exchange of fermions a with the indices (m, n) permutated. 



b) Interaction scattering ^{^1), as defined in equation (3-9), for fermion interaction between to 
composite bosons in the absence of fermion exchange, c) Exchange interaction scattering 1): 



as defined in equation (3-10), when the fermion exchange takes place after the interaction, so 
that these interactions are between the "in" cobosons d) Exchange interaction scattering 

C°"*(mi)' ^ defined in equation (3-10), when the fermion exchange occurs before the interaction, 
so that these interactions are between the "out" cobosons {m,n). 
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the "in" and "out" exchange interaction scatterings shown in Figs. |2](c) and through a 
sequence of direct interaction scattering and dimensionless Pauh scattering, according to 

rrji) = T.^{ii)KVi) (3-10) 

In these exchange interaction scatterings, the "in" and "out" cobosons are made with dif- 
ferent pairs, the fermion interactions taking place between the "in" cobosons in ^*"(mi) and 
between the "out" cobosons in ^°"*(mi)- These exchange interaction scatterings are hnked 
to the Pauh scattering through^ 

so that they are equal for energy conserving processes. Since both, interaction processes 
and fermion exchanges conserve momenta, all these scatterings differ from zero only for 

Qm ~l~ Qn Qi ~l~ Qj- 

In the following, we will see that the physically relevant combinations of energylike scat- 
terings appear to be 

iiii) = ^CA)-ru) (3-12) 

that we are going to symmeterize as 

Ln;j = liili) + im^n), (3-13) 

in order to have Cmn;ij = Cnm;ij = imn-ji- III the same way, we introduce the symmetrized 
Pauli scatterings \mn;ij^ equal to [A(J^ j) + {m ^ '^)]/2- 

B. Many-body effects 

To derive many-body effects with identical cobosons, it is convenient to iterate the four 



commutation relations (3-5 )-(3-9 ) on which the composite boson many-body theory is based. 
This leads td™ 

[5™, Ef] = iV5f -i(5^„ - D^,) - N{N - 1)5^-2 \il :)Bi (3-14) 

= ivi?r-^ y: {hi o) + (^ - n)}Bi (3-15) 
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for many-body effects dealing with fermion exchanges and 

[H, Bl^] = NBl^-\E^Bl + V^) + ^^ Bl^-' ^1 ^BIbI (3-16) 

= NBl'^-'Y.ai'DBlBi (3-17) 

for many-body effects deahng with fermion interactions. The prefactors in front of the 
A and ^ terms come from the number of ways among to choose the cobosons in state o 
involved in these scatterings. This gives when only one coboson in state o is involved and 



A^(A^ — l)/2 when two cobosons in state o are involved, the additional 2 in (3-14) coming 
from the fact that these cobosons can exchange either their fermion (3 or their fermion a, 
the (a (3) exchange being identical to a (i j) exchange, as seen in Fig. |2]^a). 

In typical problems dealing with cobosons, most of them are in the same state. Scalar 
products of such A^ coboson states are easy to expancP^ in terms of 

(0|i?f 5f |0) = N\F^. (3-18) 

The normalization factor F/v reduces to 1 for the case of elementary bosons, but it is not 
a number of the order of 1 for cobosons due to the Pauli exclusion principle between the 
fermionic components of the particles. Indeed, from the recursion relation obeyed by F^, 
it has been shownl^S that, for large samples, is exponentially small (F^r ~ e~-^^), where 



T] is the dimensionless parameter associated with density defined by (1-4). However, in 
calculations of physical quantities, F^ always appears in ratios F^_n/ F]^ which, for n small, 
is equal to 1 within corrections of the order of rj 

Fn-i 



l + 0(r/). (3-19) 



C. Coherent superposition of cobosons 

Besides many-body states with a given number of particles in a specific eigenstate, we 
wish also to consider coherent superpositions of cobosons, i.e., states of the form 

B^ = a'Bl, + a"Bl„. (3-20) 

While calculations for N = 2 are easy to perform by simply expanding {a'B\i + a" B\„Y, 
calculations for large A^ are more tricky. For those, we are forced to keep the coherent state 



i?"!" as an entity, and to produce equations similar to (3-14)-(3-17) with B^ replaced by B"^ 
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1. Fermion exchanges 

To get the fermion exchanges in a convenient way, let us first introduce the operator D^, 



through the commutator [Bm,B^ = Sm — where 6m is defined in (2-7). Dm reads in 
terms of the deviation-from-boson operators as Dm = a'Dmo' + 0'"Dmo"- From it, we then 
construct the A^n's through the commutator [Dm, = 2 ^mnB\- They read in terms of 
the Pauh scatterings as 

Amn = «'^A(J^ °,) + a"^A(^ °„) + a'a" A(J^ °, ) + A(J^ °„) . (3-21) 

Iteration of these two commutators allows us to write [Dm-, -B^^] and [Bm-, B'^^] in terms of 
this A^^.We find 

[Z}„,i?t^] = 2NB^''-'Y.>^mnBi (3-22) 

[Bm, 5^^] = NB^''-\5m - Dm) - N {N - l)fitA^-2 ^ XmnBi (3-23) 

From these two commutators, it is also possible to show that [5, 5^] = 1 — D while [D, B''] = 
2 J2 ^nB\ in which the prefactor A^ reads in terms of the Pauli scatterings as 

A„ = a'|af A(^,°;) + a"|a'f A(^„ f,,) 
+ci"\af {K, i') + XO °;,)} + a'la'f {X^,, i„) + A(^„ ii)] . (3-24) 

Similar commutators for A^, obtained by iteration, read in terms of this A„ as 

= 2iV5t^-i 5] A„5t (3-25) 
[B, B^^] = NB^^-\l -D)- N{N - l)B^^-^ J2 ^"^l (3-26) 



These two commutation relations, along with equations ( 3-22|3-23 ), are the equivalent of 



equations (3-14) and (3-15) for coherent superpositions of cobosons. 



2. Fermion interactions 

To perform calculations for of these coherent superpositions in a convenient way, 
we also need similar commutators for the interaction part. They are obtained through 
[H, B^ = EoB^ + , in which we have set = a'V^, + a" bI„ while 

B^ = {a'Eo'Bl + a"Eo"Bl„)/Eo (3-27) 
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So that B^^ reduces to for Eo' ~ Eo" ~ Eo- We are then led to define through 
B''] = J2 ^mnBl^Bj^. This scalar is found to be 

Un = o') + ci"\{l o") + a a" "J) + ^1 1")] , (3-28) 

so that the iteration of these commutators leads to 

[\/t,5tA^] = NB^''-^Y.^„,r.BlBi (3-29) 
[if,StA^] = NB^^-\E^B^ + 0) + ^^^~^^ B^^~'Y.U^BlBl (3-30) 

5. Normalization factors 

As for single composite bosons, the norm of the state made of identical coherent super- 
positions of cobosons is going to play a key role in the many-body physics of these systems. 
This leads us to introduce, just as we defined for single cobosons, the normalization 
factor Gn, determined by 

(0|5^5t^|0) = N\Gn. (3-31) 

Gat is expected to be exponentially small due to the many exchanges which take place 
between the fermions of the coherent bosons. As for F/v, we reach Gjv through the 



recursion relation it obeys. From (3-26) and the fact that -D|0) = 0, we find 

{Q\B^B^^\Q) = {Q\B^-^[B,B^^]\Q) 

= N{N - 1)\Gn-i - N{N - 1) E A„(0|fi^-^fii5t^-2|0) (3-32) 



We now use ( |3-23[ ), to get {0\B^-^Bl This leads to 

Gn = Giv-i - (A^ - 1)A(^)G^„2 + (A^ - 1)(A^ - 2)X'^^^Gn^3 + • • • , (3-33) 
where the precise value of the first order term A*^^-* is 

A^^) = E ^ns: = iar A(°; °;) + wrxc S) + 2\dar [Kt: °:) + a(°: °;,)} (3-34) 



Equations (3-33) and (3-34) show that like Fjv, the Gn correction to the bare elementary 



boson normalization factor A^!, although not of the order of 1, is such that 

1 + A^a(^) = 1 + Oil]) (3-35) 



Gn~i 
G 



N 
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D. Mixture of cobosons 



In this paper, we are also going to consider a mixture of cobosons of the form Bl^'^Bl^^\0) . 
Like F/v for many-body effects between identical cobosons, Gn-^^n2j defined by 



Ni,N2, 



(3-36) 



is a key factor for many-body effects between Ni cobosons in state Oi and N2 cobosons 
in state 02. This factor, which would be exactly 1 for elementary bosons, differs from 1 
due to the many fermion exchanges which exist not only among the A^^i cobosons in state 
Ol and among the N2 cobosons in state 02, but also between the cobosons in state oi and 
the cobosons in state 02- They make Gni,N2 not of the order of 1 but exponentially small. 
As with F/v, we can determine Gni,N2 through the recursion relation it fufills. To get this 



recursion relation, we use equation (3-14) to rewrite Bo^Bl^^. This leads to 



(3-37) 



The 1 in the second term in the bracket readily gives Ni [{Ni — 1)\N2\Gn-^-i^n,j]. The first 



term in the bracket is calculated using (3-14) for Ba^Bl^'^\0), while the third term is calcu- 



lated using (3-15) for Doid-B^j^^ |0). This allows us to write the norm of Bl^'^Bl^'^\0) given 



in equation (3-36) as 



N2\N2\Gn,,N2 = iVi[(iVi - 1)\N2\Gn,-i,N2] - N,{N, - l)Au - N2{N2 - 1)^22 - NiN2A,2, 

(3-38) 

where An, A22, and yli2 contain one Pauli scattering explicitly, so that they are related to 



fermion exchanges. Equation (3-38) is shown in Fig. the factors {Ni,N2) coming from the 



number of ways we can choose the cobosons involved in fermion exchanges with the coboson 
Ol on the left. The leading term for A'^i(iVi — l)Aii, shown in Fig. Ill reduces to 



(iVi-l)A(°^ :i)N^iN,-l) 



{Ni-2)\N2\G 



Ni-2,N2 



(3-39) 



Similarly, the leading term for A'^iA^2^i2, shown in Fig. |5| is 



No 



02 01 ^ 
Ol 02 1 



NiN2 



(N^-l)\{N2-l)\GN-,-l,N2-^ 



(3-40) 
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FIG. 3: Shiva diagram expansion for the normahzation factor {0\B^^ B^^ Bl^^ Bl^^\0) appearing 



in equation ( 3-38 ) . 
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FIG. 4: Shiva diagram expansion for the first exchange contribution iVi(A'^i — to 



(0|S^2 5^i5tiVi^tAf2|o), as given in equation (3-39). 
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FIG. 5: Same as Figure 4 for the third exchange contribution A^iA''2^i2, as given in equation (3-40 1 . 



They contain fermion exchanges between two cobosons. On the opposite, the leading term 
of N2{N2 — 1)^22 niust have fermion exchange between 3 cobosons since X{°^^ 
0\ 7^ 02. Therefore, we end with 



°2) = for 

02 02/ " 



Gni,N2 — Gn^-I^N2 - {Ni - °o\)Gn^-2,N2 - N2 °o\) + \{°o\ 02)) Gn^-I^N2-1 + 



(3-41) 



(3-42) 



This shows that, in the same way that F/v-iZ-F/v ~ 1 + 0(77), we have 

^^^ = l + (9(,,) + (^M, 

^Ni,N2 

SO that, even if Gni,N2 is exponentially small due to the large number of fermion exchanges, 
the effect of these exchanges on Gni,N2 ratios are negligible at lowest order in density. 

We now have all the tools to tackle many-body effects with a large number of identical 
cobosons, a large number of coherent cobosons and a mixture of two large numbers of 
cobosons. However, since the calculations for large A^'s are obviously quite technical, in this 
paper we have chosen to start with = 2, as most of the many-body physics can usually 
be understood from this limit. We are also going to first perform these calculations for 
elementary bosons, in order to enlighten the differences between elementary and composite 
particles. 
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IV. TWO ELEMENTARY BOSONS 
A. Single state 

Let us start with a state having two identical elementary bosons, |0o) = Bl'^\0). Since 



Bi\(f)Q) = 26ioBl\0) due to equation (2-4), we find that (0o|0o) = 2 while 



(0o|iyo|0o) = Y.E^{{h\B^}{BSo)} = 4^o. (4-1) 

(Mvi^o) = lj:^c:ni){^o\BiBiB,BM = 2^::). (4-2) 

This readily gives the Hamiltonian mean value in this two-elementary boson state as 

" \H\^o) 



^01 



(0o|0o) 



2Eo + aro) = ^Eo + U,oo, (4-3) 



where ^mn;ij is defined in terms of ^(^ }) as in (3-13). 



B. Fragmented state 

We now consider |0i2) = BI_^BU0) for elementary bosons in states Oi 7^ 02. From 
Bi(j)12), we readily find (0i2|0i2) = 1, while (0i2|-f^o|0i2) = -E'oi + -£'02- Iii the same way, 
BjBi\(l)i2) = {Sjo^Sioi + Sjo^6io,,)\0) leads to 

{^i2\v\^i2) = a:roi)+a:i:i), (4-4) 

so that we end up with 

(0121012) 

for oi ~ 02 ~ o. 



C. Coherent superposition 

The third state of interest is the coherent superposition of states, |0) = 5t^|0), with B^ 
defined as in (3-20), with bI, replaced by -Bo'^- The prefactors {a', a") are chosen such that 
|a'p -|- |a"p = 1 in order for the state to be normalized. We also impose that the {o',o") 
states differ at least through their center of momenta {Qo' 7^ Qo") in order for scatterings 
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like ^(°'' °I) to cancel due to momentum conservation. From equation (2-6) taken for N = 2, 
we readily find = 2 while (0|ifo|0) = 4(|a'pi?o' -|- \a"\'^Eo"). In the same way, from 

BjBi\(j)) = 2{a'6io' + a"6io"){a'6jo' + a"6jo")\0) and the fact that the nonzero ^ reduce for 
Qo' 7^ Qo" to the diagonal terms ^(* *•), with (i, j) = (o' or o"), or to the cross term ^(°I, °r), 
we find that 

(01^10) = 2\amo' o') + 2|a're(o" o") + 4|a'a'f [e(°;' °;') + ^o" o')] ■ (4-6) 

Consequently, the Hamiltonian mean value for two elementary bosons in a coherent state 
reads as 

((t)\H\(t)) /|2 771 I I "i2t7 \ I I /|4^^ I I //|4^^ I Ai I n\2t 

- Z[\a I r^o' + P I ^o") + P I ^o'o';o'o' + P I ';o"o";o"o" "h 4|a I .^o'o";o'o" 



(010) 

^ 2Eo + (1 + 2|aVf )U;oo (4-7) 
for o' ~ o" ~ o. 

From equations ( 4-3|4-5|4^ ), we thus see that, since the interaction scattering ^oo;oo must 



be positive to avoid a density collapse, the energy of two elementary bosons is minimum 
when the particles are in a single eigenstate, and not a mixture or a coherent superposition 
of two eigenstates. 

V. TWO COMPOSITE BOSONS 

The calculations for composite bosons will need to explictly treat the exchange between 
the constituent fermions, through the deviation-from-boson operator Dmi and the Pauli 
scattering A (J: 



(n 



A. Single state 



For |0o) = -Bq^IO), we now have, due to equation (3-14), 

Sn.|0o) = 25^oi?t|0) - 2 ^ AC^ ^)St |0), (5-1) 
since -0^1 0) = 0. This leads to 

(0o|0o)=2-2A(ro) (5-2) 
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with A(°°) = (337r/2)(aB/L)3 for 3D cobosonpl with zero center-of-mass momentum and 
relative motion wavefunction (r|z/o) = e~''/"^/y^7ra^; so that this Pauh scattering goes to 
zero when the sample size L increases. 



Using equation (3-16) and the fact that Vg |0) = 0, we find 

H\ct>o) = 2E^Bl'\Q) -Y.ai:)BlBi\Q) (5-3) 

From equations ( 3-6[3-7 ), it is easy to show that the scalar product of two coboson states 
reads a^ 

{Q\BmBr.B\B]\{)) = {5ra^5n, ' Km ))} + {m ^ u} (5-4) 



So that, from the two above equations, we end up with 

(0o|i^i0o) _ ^ e(°g)-r(g°) ^ + 



1 ^oo\oo 



(0o|0o) l-A(n) 
since \oo;oo goes to as {aB/LY when the sample size increases 



2£'o ~l~ 'Coo;c 



(5-5) 



B. Fragmented state 



We now turn to the fragmented state |0i2) = Bl^B\^\Q) for (oi, 02) with different center of 
mass momenta in order for the Pauli and interaction scatterings to cancel due to momentum 
conservation if the number of Oi cobosons in the "in" and "out" states are different. Using 



eq. (5-4), its norm is found to be 



^121012) = i-KlZ)-K\:i] 



1-2A, 



oi02;oi02 ) 



(5-6) 
(5-7) 



while from equations (3-8)-(3-9), the Hamiltonian mean value in this fragmented state ap- 
pears as 



^i2|g|0i2) , e(gn^)-r(gg^) + e(sg)-r(g^g) 

(0121012) °' 1-2A 

0102,0102 



(5- 



So that, if we use the fully symmetrized scattering defined in (3-13) and take a large sample 



volume in order for the A term to give a negligible contribution, this Hamiltonian mean 
value reduces, for (01,02) ~ o, to 



''12 



\H\ 



^12 



) 



(0121012) 



Eo^ + + 2 



^01 c 



02 



A, 



(5-9) 
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C. Coherent superposition 



The third state of interest is the coherent superposition |0) = B'^'^\0), where B"^ = a'B\, + 
a"Bl„, with |a'p + |a"p = 1 and Qo' 7^ Qo" in order for the Pauh and interaction scatterings 
to again cancel due to momentum conservation if the number of o' cobosons in the "in" and 
"out" states are different. As a direct consequence, and {(f)\H\(j)) can only have terms 



I /I4 , .. 

m a , a 



bI,bI„, readily leads to the norm of the coherent state given by 



"1^ and |a'a"p. Using (5-4), the bare expansion of -B^^ in terms of bI?, B^?, and 



/ 1 4 r I I " 1 4 r I I ' " 1 2 r 

CL L/qI -\- \(l L/qII -\- \(1 d LiqI q: 



(5-10) 



with Lq/ = 2 — 2A(°I °J) and similarly for Lg", while Lq/^o// = 1 — A(o" °") — A(°» °"). In terms 



of the symmetrized Pauli scatterings defined as in (3-13), this norm reduces to 



2 — 2|o'| Xq'o'-o'o' — 2|ci"| \o" o" -o" o" — 8|ct'a"| \o' o" o' o" • -'--'-) 



In the same way, (3-8) and (3-9) allow us to write 



where An' = 2En' 



{(t>\H\(l)) = \a'\^Ao' + \a"\''Ao>' +4\a'a"\'Ao',o", 
2 - 2A(°; + 2^"(°; °;) and similarly for Ao", while 



(5-12) 



Ao',o" = {Eo' + Eo") 1 - A(°, °, ) - A(°„ °, ) + ^(°, °, ) + ^(°„ °, ). 

For o' ^ o" ~ o in a large sample volume, the expectation value of the Hamiltonian then 
reduces to 



\H\ 



2Eo + 



1 + 2\a'a" 


J<,oo;oo 




2 


a'a" 


/ ^^oo;oo 



2Eo + (l + 2|a'a'f ) Lo;oo (5-13) 



D. Effective scattering for bosonized particles 



When compared to similar results for elementary bosons, namely equations ( 4-3[4-5[4-7 ), 
the Hamiltonian mean values in the three two-coboson states, obtained above, lead us to 
identify the diagonal effective scattering for elementary bosons ^oo-oo with the physically 



relevant combination of energylike scatterings defined in (3-12), namely 



C,oo;oo — C,oo;oo — ^{o o) ^ \o oJ- 



(5-14) 
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This effective sacttering contains a direct contribution as well as an exchange contribution 
which is symmetrical with respect to the "in" and "out" processes, as physically reasonable 



for f "(o o) = r"*(o o) due to equation (3-11 ). 



For excitons or H atom^, the diagonal direct scattering °) reduces to 0, the repulsion 
between fermions a or fermions /3 being as large as the attraction between On the 

opposite, the diagonal exchange scattering of these cobosons differ from zero: in 3D, it reads 
= -(267r/3)i?o(aB/L)3, where i?o = [le^ll^h^e^ = jlas is the coboson Rydberg 
energy.!^ 

VI. MANY ELEMENTARY BOSONS 

We now turn to states with a large number of bosons and first consider that these bosons 
are elementary bosons. 

A. Single state 

Let us start with the state having its bosons in a single state, |0o) = -Bj^|0). We 
physically expect the interaction term of the Hamiltonian expectation value to depend on 
the boson number through N{N — l)/2 which corresponds to the number of interacting 
boson pairs (o, o) we can form out of bosons in state o. This leads us to expect that 



the Hamiltonian expectation value obtained for = 2, as given in equation (4-3), must 
transform for general as 

(^o^l^o) =NE+ (6-1) 
For scatterings ^oo;oo in {(^b/LY, this will induce a correction to the bare energy NEo of the 



order of Nr], where r] is the dimensionless parameter associated to density given in (1-4). 
Let us now recover this physically expected result. 



Equation (2-4) readily gives the well known normalization factor for elementary bosons, 
namely 

(0o|0o) = (0|B„^-iB,sr |0) = Ar(0|5„^-^Sf -^|0) = AT! (6-2) 



Using the same ( 2-4[ ) on the one-body part of the Hamiltonian leads to 



^o|^o|0o) = E Ei{{^o\Bj} {b,\^o)} = N^Eo{N - 1)! = ArEo(0o|0o). (6-3) 
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while on the two-body part, we find 

{Un^) = \Y.'i{Vi)[{MBlBi}[B,Bl^,)} = lN\N-in{::){N-2)\ (6-4) 



This just proves that the guess of equation (6-1) is correct. 



B. Fragmented state 

We now consider the fragmented state \(l>i2) = Bl^'^Bl^'^\0). The same physical under- 
standing leads us to expect the interaction term between elementary bosons (01,02) of the 



Hamiltonian expectation value, as given in (4-5), to appear with a N1N2 prefactor. Since 



we now have many bosons Oi and many bosons 02, we should also have an interaction term 
between bosons oi with a prefactor Ni{Ni — l)/2 and an interaction term between bosons 



02 with a prefactor N2{N2 — l)/2. From(4-5), we are thus led to guess 

(012|i^|012) ^ 

+2 A^l A^2'Coi02 ;oi02 



(0121012) 



NF ^NF I ^1(^1 -1) .- I N2{N2-l) , 

-''l-'-'oi ~r -'»2-'-'02 ~r d Soioi;oioi ~r „ S0202;0202 



(6-5) 



For oi ~ 02 ~ o, the Hamiltonian expectation value should thus be given by 



n2 



1^1012) 



11 



(0121012) V 2 

for Ni + N2 = N. Let us now show this result explicitly. 



+ NiN2 U;c 



(6-6) 



Equation (2-4) leads to 



5,1012) = (N,6,^M^^-'Bl^'+N2S.oM^^Bl^'-') |0) 



(6-7) 



so that, for oi 7^ 02, we find 5^i|0i2) = Ni\Bl^^\0). This readily gives (0i2|0i2) = A^i!A^2!, 
while 

(0i2|i^o|0i2) = Y.M^i2\BjBSi2) = E,^N^{N, - 1)\N2\ + E,^Nl{N2 - l)!iVi! (6-8) 



If we now use equation (6-7) to calculate -Bj-Bj|0i2), the scalar product |(0i2|-B,^-Bi| 



|i?ji?i|0i2)| which appears in (0i2|^|0i2), leads to 



(6-9) 
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in which we have set 



= a:\ :\)[Ni{N^ - i)f{N, - 2)m 



and similarly for W021 while 



The expected result (6-5) then readily follows from equations (6-8)-(6-9). 



(6-10) 



(6-11) 



C. Coherent superposition 



The third state of interest is the coherent superposition of states |0) = B'^^lO) with 
B'' = a' B\, + a" B\„. Through similar physical arguments, we expect to have the interaction 



term in the energy for N = 2 cobosons, as given in (4-7), to appear with a prefactor 



A^(A^ — l)/2. To show it explicitly, we use equation (2-6) to get 



{0\B^-\a'*Bo' + a"*Bo")B^^\0) = N[\a'\^ + |a'f ](0|S^-^5^-^|0). (6-12) 



Its iteration, for |a'p-|- |a"p = 1, gives = N\ as for a condensate made of a single state. 



In the same way, the norm of Bi\(j)) deduced from equation (2-6), leads to 



Ho\(j)) = N{\a'\'Eo^ + \a"\'Eo")N\ 



(6-13) 



while the scalar product of the states {(f)\Bl^Bl and BiBj\(f)) leads, for Qo' 7^ Qo", to 

= \[N{N - l)nN - 2)! {\a'mi, i,) + |a'f e(o" o") + 2|a'a'f (e(°;' ^j') + °;'))} , 

(6-14) 



which also reads 



{^\V\4>) = ^^^o [W\%'o';o'o' + |a'f eo"o";o"o" +4|a'a'f eoV';o'o"]A^! 



So that, for o' ~ o" ~ o', we end with the expected result, namely 



(010) 



NE. 



(6-15) 
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N 



FIG. 6: Shiva diagram for a process having one fermion interaction between cobosons o, as given 



by the second term of equation (7-1 ). 

VII. MANY COMPOSITE BOSONS 

A. Single state 

The Hamiltonian mean value for the pure state |0o) = Bl^\0) has already been calculated 
in a previous worlpSESl. it has a naive contribution NEg. It also has a set of density 
dependent corrections in rj'^ with > 1, in contrast with elementary bosons which only 
have a n = 1 term. This set of density terms comes from fermion exchanges between the 
N cobosons: Since the Hamiltonian expectation value (H) only contains one interaction 
by construction, the density terms for n > 2 can only come from the fermion exchanges 
between 3 or more cobosons induced by the Pauli exclusion principle. They are nicely 
visualized by Shiva diagrams with n + 1 cobosons and one interaction process between any 
two of these coboson lines (see Fig. [6]). Let us here repeat the main steps of this calculation 
for completeness - and also because the ones for \(pi2) and |0) are conceptually similar, while 
far more complex. 

In order to calculate {(f)Q\H\(j)Q), we use one of the two key equations for coboson many 



body effects, namely (3-16), to find 



((/.o|i/|0o) = iVi?o(0o|0o) + ^^^^ ^^ E^rmimB^Bl'^-'BlBilO). (7-1) 

The second term of this equation is shown in Fig. [6] In it, appears the scalar product of 
N coboson states with two cobosons different from o on the right (see Fig. [T]). This scalar 



product is calculated using another key equation for many-body effects, namely (3-14). This 
equation leads to 

(OlSr^S^S.^nO) = {0\B^-'B,^ [NSr.oBl''-' - NiN - 1)5^"' E ^(i; o)^;] |0). 

(7-2) 
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FIG. 7: a) Shiva diagram for the scalar product of N coboson states appearing in equation (7-2), 
with N cobosons o on the left and (N — 2) cobosons o plus two cobosons (m, n) on the right. The 
standard to calculate this scalar product is to isolate {N — 2), (A^ — 3), . . . cobosons o not 

involved in fermion exchanges with (m,n). The prefactors come from the number of ways to 
choose the coboson o having fermion exchanges with {m,n). These possible exchanges are shown 
in (b) and (c). 
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We then use the same commutation (3-14) to pass Bm over Bl^ ^ and Bl^ ^. Since 



■^(m }) — Hi m)' this allows US to expand the above scalar product as (see Fig. |7| 
{0\B!^BlBiBl^-'\0) = 

N{N - 1) [6onSom - K ™)] (N - 2)\Fn-2 

+NiN - l)(iV - 2)2 [-5„„A(° - 6omXi: ™) + A3 + A'3] (N - 3)!F;v-2 + 



where A3 and A3 are three-leg scatterings shown in Fig. [Tj^c). 



(7-3) 



A3 = A 



o n 
yo m J 



A'3 = A 



^ o 



o o 
yo m J 



in which the cobosons (m, n) exchange their fermions with one coboson o to produce three 
cobosons o. This expansion actually follows the standard procedurd^ to calculate scalar 
products, namely, we first isolate terms in 

{0\B^-''Bl''-''\0) = {N- P)\Fn-p, 

with P >2. We then connect the remaining cobosons in all possible ways while enforcing the 
cobosons o to be "never alone" as in Fig. [7]^b),(c). The N prefactors in Fig. [Tj^a) correspond 
to the number of ways we can choose the cobosons o among on the left and among 
(A^ — 2) on the right. In the case of the first term of this figure, we just have to choose the 
two cobosons o on the left; this is why this first term appears with a prefactor A^(A^ — 1). 
In the second term, we also have to choose a coboson o on the right and a third coboson o 
on the left. This is why this term appears with a prefactor [N{N — 1){N — 2)](A^ — 2). And 
so on ... 



By inserting equation ( |7-3[ ) into ( |7-1[ ), we readily find 



NEo + 



N{N 



e(oo)-e 



o o) 



+ ... 



so that, since Fisf-2/FN = {Fn_2/Pn-i){Pn-i/Pn) = 1 + 0{ri), we end for large A^ with 



(7-4) 



\H\<Po) 



) 



A^ 



(7-5) 
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This shows that the first correction to the bare energy Eo is of order t] = N{aB/L)'^, 
because, as shown in Section V.D, C,oo]oo is of the order of {as/LY. Note that, as for 
Coulomb interaction = 0, the effective scattering ^oo;oo reduces to —C"'{oo) which for 

3D excitons is equal to —(267T/3)Ro{aB/LY. Consequently, this Hamiltonian mean value is 
just the energy of N electron-hole pairs obtained by Keldysh and Kozlo'v'^, using a completly 
different approach in which these pairs are not treated as coboson entities, as we do here. 



B. Coherent superposition 



In order to calculate the mean value of the Hamiltonian for the N coherent cobosons 



-B^^lO), we use the results derived in Section IIIC From (3-30), this mean value 
appears as 

... . N(N - n 

(7-6) 



(7-7) 



where ^mn is given in eq(3-28). The first term of (7-6) reduces to NEo{(f)\(l)) for Eo' ~ Eo^ 



Eo since B'' ^ B'^ . To calculate A, we use equation (3-26) for {{)\B^ B}^. This leads to 



In the first term, we again use (3-23) for {0\B^ ^i?^, while we use the commutator [Bp, Bj^] 



given in (3-6) to calculate the second term of A. This leads to 
A = iV(iV - l)(iV - 2)\GN-2j2UniS*j: - - 



...Gn-3 + 



(7-9) 



If we now use the definitions of Sm, imn and A^n given in equations (2-7), (3-21 ), and (3-28) 



we find that for Q^' 7^ Qo", the sum in (7-9) reduces to 



with ^ = ^ — as defined in equation (3-12). By collecting all the terms, we end with 
N{N-1)Gn-2 



NEo + 



G 



N 



{k'I^Oo';o'o' + W\'^io"o";o"o" + 4:\ct' d" \'^io'o" ^o'o"^ + 



(7-10) 



within corrections of the order of {Eo' — Eo"). For o' k, o" ^ o and N large, the above result 
reduces to 

(7-11) 



mil 

(010) 



N[Eo + + 2|a'a"nU;oo + 0\v% 
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C. Fragmented state 

The last N-coboson state we must consider |0i2) = Bl^'^Bl^^\0) has two large numbers 
of different cobosons in state Oi and 02. The very many exchanges which exist within the 
Ol population, within the 02 population, and between the oi and 02 populations make this 
many-body calculation quite tricky. This is why we have kept it for the end, the previous 
ones having the role of useful exercises. In view of the above results, we can however guess 
that the Hamiltonian mean value for composite bosons in this fragmented state should read 



as the one for elementary bosons, namely (6-5), with, according to (5-14), ^oioi;oioi and 
^oi02;oi02 replaced by ^oioi;oioi and ^0102,0102: respectively. Let us show this nicely simple 
explicitly. 



We first use (3-16) in the coboson many-body effect section, to replace HB^^^ in if |0i2) 



This leads to 

H\<t>.2) = I^B^^^H + N,E^,Bl^^ + N,Bl^^-WX + ^^^^^ ~ i^f-^ E IDBIbI^ 5^10). 

(7-12) 



We again use (3-16) for HBl^^\0) and (3-17) for V^Bl^^lO). This allows to split (0i2|if|0i2) 



into 

NUN, - 1) NoiNo - 1) 

(012|if|012) = (NiEo, + N2Eo2){h2\h2) + ^An + ^A22 + iViiV2Ai2. 

(7-13) 

All, shown in Fig. |8j describes the interactions of two among A^i cobosons in state Oi, the 
other N2 cobosons 02 possibly having fermion exchanges with the cobosons Oi. The precise 
value of All is 

An = E mBSK'BlBiBl^^-'Bl^^mH ol), (7-14) 
and similarly for A22, while A 12, shown in Fig. |9} results from the interaction of one coboson 
Ol with one coboson 02. It reads 



A12 = E {OK^'B^^^BlBiBl^^-'Bl^^-'mi 



n 02 ^ 
m 01 y 



(7-15) 



The reader who is knowledgable about Shiva diagram^^MI] -^jj] immediately see from 



Fig. 8 that, for Q^, ^ Q^^, i.e. for ^("^ °;) = = f "("^ °i), we must have 



Aii~iVi(iVi-l) 



VOl Ol, 



-c 



{N^-2)\N2\Gn,-2. 



N2, 



(7-16) 
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FIG. 8: Shiva diagrams for An defined in equation (7-14), in which two cobosons oi among Ni 



have an interaction scattering, the other cobosons 02 just possibly exchanging their fermions 
with the cobosons oi . The prefactors are the number of ways to choose the cobosons oi on the 
left. 

and similarly for A22- Here Gjy-^^iy^ is defined as Fjy, through the norm of the Bl^^Bl^^\0) 



state (see equation (3-36)). 



In the same way, the Shiva diagram of Fig. M readily leads to, for Qo^ 7^ Q 



,02 1 



A 



12 



NiN, [e(^j°j)+e(s°?)-r(o?o?)-r(o^o?) {n,-i)\{n,-i)\Gm,^,,m2^,. (7-17) 



By collecting all these terms, we end with 



'12 



'12 



iVi(iVi-l)^ 



Gn,-2,N2 , N2{N2-l)^ 



G 



'12 



1012) 



— -^lEoi + N2E02 "I" r) ^OlOi;OlOl ~^ r, ^0202,0202 x-y 



Ni,N2-2 



Gni-1,N2-1 

02\Ol02 

'^Ni,N2 



Ni,N2 

(7-18) 



Like F/v, the factor G'7Vi,Af2 5 which comes from the many fermion exchanges which take place 
between the A^^i cobosons in state Oi and the N2 cobosons in state 02, is exponentially small. 
However the ratios of Gni,N2^s are nearly 1 at lowest order in density. Enforcing A''i-|-A^2 = N 
and oi ~ 02 ~ o, we ultimatly find the expected result, namely 



'12 



'12 



(7-19) 



(0121012) - - \ 2 

Let us end this section by returning to the interaction parts An, A22 and A12 of the 



Hamiltonian mean value defined in equations (7-14) and (7-15) and by calculating them, 
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FIG. 9: Same as Fig. 8, for A12 defined in equation (7-15), the interaction taking place between 



one coboson 01 among A^i and one coboson 02 among The prefactors are the number of 
ways to choose the cobosons (oi, 02) on the left. 

not through Shiva diagrams, but in a pedestrian way using the commutators appropriate to 



many body effects, namely equations (3-14) and (3-15). The problem is to get the scalar 



products of coboson states which appear in these A's. This always is the tricky part of any 
calculation involving cobosons. The ones of interest here are 



(7-20) 



for P = (0, ±1). This calculation is done along a line similar to the one we have used when 
we only had one type of cobosons, namely we isolate the norm of states like B!,,'B^r\0) 
with {N[,N2) decreasing from {Ni — 1 — P, N2 — I + P). To do it, we first use equation 



(|3-14|) to rewrite B^^^Bj^. This leads to four terms. We then use (|3-14|) to get {0\B^^^B}^ 



and equation (3-15) to get {0\B^^Dom- This allows to split S^^^ as 



= NiTP + N2TP - Ni{Ni - 1)U^P - N2{N2 - l)a 



(P) 



riP) 
22 



12 ) 



(7-21) 
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where the contributions Ti^'' and T2^^ are somewhat direct since they read 



n 



where we gave set |Vn,p) = 5t5t/i-i-^5t/2-i+P|o) 

(p) (p) (p) 

The three other terms f/ji , U22 and U12 contain one Pauh scattering exphcitly which 
describes the fermion exchanges between two cobosons in state oi, two cobosons in state 
02, and one coboson in state oi with one coboson in state 02, as understood from their 
prefactors. These f/*^^-* terms precisely are 



(7-22) 
(7-23) 



tAP) 

^11 



MP) 



MP) 
'^12 



Ol t 
oi m 



EA( 

i 
i 



02 1 

02 m 



B^r'K'B.\i^n,p) 



B^r'Kmn,p) 



(7-24) 
(7-25) 
(7-26) 



Their leading contributions in fermion exchanges are obtained by passing Bi over Bj^ in \ipn,p) 



through the commutator (3-6). The trivial term corresponds to taking i = n while the two 
other terms generate additional exchanges between i, or n, and the other cobosons in states 
Ol and 02. In the case of u[i \ the remaining matrix element is just {Ni — 2)\N2\Gn-^-2,N2 
for P = 1 while it contains additional Pauli scatterings for P 7^ 1. In the same way, the 
remaining matrix element of UiP is just {Ni — 1)!(A^2 — l)'G'iVi-i,Af2-i for P = 0, while 
for other P's, this matrix element contains additional Pauli scatterings. Consequently, the 
contributions to the f/^^^'s with only one Pauli scattering reduce to 



^11 



MO) 
^12 



A(°; 1){N,-2)\N2\G,,,^2,N2 
KrJ + Kll)] (iVi-l)!(iV2 



(7-27) 



Ni-1,N2-1 



and similarly, U22 obtained from by changing 1 into 2. 

AP) 



(1) 



Let us now turn to the Ti term defined in (7-22). For P = 1, we pass B^^ over 



Ni-l 



using the many-body commutator (3-14). The contribution without exchange, which is the 



dominant one at small density, reads 

) ~ (iVi - l)6o,n,6oUNi - 2)!iV2!G;vi-2 N2- 
If for P = 0, we do the same but with P^^, we find 

~ N26o,mSo2n{N, - 1)\{N2 - 1)!G^,-1 iV.-l. 



(7-28) 



(7-29) 
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On the opposite, for P = —1, additional exchange processes are necessary to transform the 
BI^^ operator on the right into 5^^^^^^ in order to get rid of the B^^~^ operator of the left, 
as necessary to generate a G'iVi-i,Af2 factor. By calculating T2 along the same line and by 
inserting all these matrix elements into An and A12 given in equations ( 7-14|7-15 ), we end 



with the Hamiltonian mean value written in equation (7-18). 



VIII. DISCUSSION 



In the previous sections, we have performed detailed calculations of the Hamiltonian mean 
value for three different types of states, namely a pure state |0o) = -Bl^|0), a coherent state 
10) = {a'Bl + a"Bl„Y\Q) and a fragmented state |0i2) = Bl^^Blf'\Q) with A^^ + A^2 = N. 
We have taken N = 2 first, since the calculations are rather trivial and most of the physics 
can already be understood from this two-body problem. As explicitly shown, the results for 
large N can be deduced from the ones for = 2 by putting N in front of the scatterings, in 
this way producing a density correction in 77 = N{aB/LY to the bare energy, as physically 
reasonable. 

We have considered elementary bosons as well as composite bosons. The calculations with 
composite bosons are done using the new many-body theory designed for them, in which 
the fermion exchanges appear explicitly through a set of Pauli scatterings Xi^l). These 
exchanges are visualized by Shiva diagrams, which are of great help to drive the algebra 
in the right direction in case of many body effects - the last section being quite convincing 
with respect to the utility of these diagrams. The Hamiltonian used when the composite 
nature of the particles is retained, is the microscopic Hamiltonian for fermions. On the 
opposite, an effective Hamiltonian is necessary when the composite bosons are replaced by 
elementary particles. In order to recover the exact composite boson results from this effective 
Hamiltonian, for the simple problem in which only one interaction scattering appears — the 
Hamiltonian mean value being first order in the interaction by construction — we must adjust 
the diagonal scattering of this effective Hamiltonian to be such that 

where is the symmetrical combination of direct and exchange interaction scatterings of 
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two composite bosons, defined as 



1 



1 r 



= 2 [ea I) + a) i)\ - 2 ^) + ^^"0 1)\ ■ (8-2) 

The term ^(^i) corresponds to the interaction scattering of the composite boson many- 
body theory which describes fermion interactions without fermion exchange, while the term 
C"(mi) = S -^(m p)^(p i ) corresponds to the exchange interaction scattering, with A(^^) 
describing fermion exchange in the absence of fermion interaction. Note that, in the case 
of diagonal scattering, we do have ^""(11) = i); physically required by the time 

reversal of the scattering. 

Within this identification of the diagonal effective scattering for bosonized particles, we 
find that the Hamiltonian mean values read the same for elementary and composite bosons 
in the three states of interest, namely 



^>o|0o) 2 



. , , _ NEo + ^(1 + 2\a'a"\')Uoo (8-4) 

(010) 2 

..^E„J^^^N,N^i^, (8-5) 



The above results are exact for elementary bosons and only approximate for composite 
bosons, being valid to lowest order in density only: indeed, in the case of elementary bosons, 
the many-body physics is induced by 2 x 2 interactions while for composite bosons, a quite 
subtle new set of many body effects arise from fermion exchanges which can exist between 
more than 2 cobosons. 



The single-state case is recovered for a' or a" = in ( 8-4 ) and for A^^i or A^2 = in ( 8-5 ) , as 
expected. We also see that, since ^oooo must be positive (otherwise the system would suffer 
a density collapse, its energy decreasing with increasing density), the minimum energy is 
obtained for macroscopic occupation of the single state |0o)- Consequently, a condensate in 
just one state is stable — a condensate fragmented into two eigenstates or in a superposition 
of different eigenstates has a macroscopically higher energy, even when the other states 
(oi, 02) or (o', o") are infinitesimally close to the ground state o. 

We end with one last comment. We have here considered the possibility of the condensate 
to differ from a macroscopically occupied single state through the study of the Hamiltonian 
mean value for the states |0o)) |0) and |0i2)- Of course, Bl^\0) is not the true ground 
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state of the system. It however corresponds to this ground state to lowest order in the 
interaction, as can be seen from the fact that the Hamiltonian mean value of this state has a 



zero-order term in scattering which is the expected one, namely NE^ (see equation (8-3)). It 
is important to stress, however, that while "lowest order in the interaction" is a well-defined 
concept in the case of elementary particles, for which the Hamiltonian reads Hq + V, such a 
concept has no clean meaning for composite bosons, since it is not possible to describe the 
interaction between cobosons as a potential, and thus define a zero-order Hamiltonian or a 
zero-order eigenstate. 
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IX. APPENDIX: NOZIERES'S ORIGINAL ARGUMENT 

For completeness, we here reproduce the original Nozieres's argument on Bose-Einstein 
condensatd^, that there is a macroscopic number of elementary bosons in a single quantum 
state, not from a difference in the free particle kinetic energy but the exchange part of the in- 
teraction energy. Nozieres calculates this interaction energy, within the Born approximation, 
namely 

for a structureless interaction Hamiltonian 



2 ^ ^ ki k2 «3 fel+K2— «3 

ki,k2,k3 



He first considers an elementary boson condensate made of a sing le sate |0o) = bI^\0), the 
elementary boson operator Bj: being such that [Bj^,Bp] = 5^^. The interaction energy is 
then found to be, in the large N limit, 

{V), ^ ^VoN'. (9-7) 

Nozieres then considers a condensate made, not of a single quantum state, but of two de- 
generate or nearly-degenerate states, \4>i2) = -Bi^^-B^^^ |0), with the same total number of 
particles A^^i + A^2 = N. In this fragmented state, the interaction energy is found to be 

{V)i2 = {V)o + VoN,N2. (9-8) 

For Vo positive, as necessary to prevent a density collapse, this readily shows that it we 
must pay a macroscopic amount of energy to break up the condensate into two parts, due 
to exchange between these two parts, as seen from the N1N2 prefactor of this additional 
energy. 

Instead of the structureless constant scattering Vq used by Nozieres, we have, in this 
paper, decided to use a scattering ^{^1), which a priori depends on the "in" and "out" 
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states. The idea is to make an easier comparison with the results obtained for composite 
bosons which read in terms of the specific combination of direct and exchange processes given 



in (3-12), namely ^(^D = ^(mi) ~ C"(mi)- This, in particular, allows us to identify the 
proper effective scattering we must use for bosonized particles. Let us however stress that 
this identification can only be done for energy conserving processes as in the case of diagonal 
processes, i.e., processes in which the "in" and "out" states are identical. Indeed, such a 
scattering ^(^ i) cannot be used in general because, as = JT)' ^^e resulting 



Hamiltonian would not be Hermitian, due to (3-11) which tells that the "in" and "out" 
scatterings are equal only when energy is conserved. This difficulty is basically linked to the 
fact that there is no way to have an effective Hamiltonian for bosonized particles which is 
Hermitian and valid for all many-body effects. By taking a real constant Vq, Nozieres hides 
this difficulty. 



